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ABSTRACT 

The extended and X-ray emitting interstellar medium of early type galaxies is often 
used as a tool to determine their total mass M and stellar orbital anisotropy (3 profiles, 
based on the hypothesis of hydrostatic equilibrium for the hot gas. Here we investi- 
gate the effects that deviations from equilibrium have on M and P estimates, by using 
simple analytical calculations and hydrodynamical simulations representative of gas 
rich galaxies. We show that the deviations of the X-ray determined Z?*^** and M*^** 
from the true values are linked by a remarkably simple relation; in particular, M is 
underestimated if /? is overestimated. Also, more radially anisotropic orbital distribu- 
tions than true are deduced in presence of gas infall velocities of the order of the local 
stellar velocity dispersion (as are likely in the central regions of galactic cooling flows). 
The results of this analysis are applied to the most thoroughly investigated bright el- 
liptical, NGC4472. First we show that /3^''* recently derived from X-rays corresponds 
to a galaxy that is unstable by radial orbit instability. Then, assuming as true (3 and 
M the optically derived values, we show that the differences /J*^"* — f3 and M^** — M 
agree with the predictions found here in the case of lack of hydrostatic equilibrium, 
which points to the latter as a possible explanation for the discrepancies. This analysis 
casts doubts on the possibility to use X-ray information to determine accurately the 
dynamical properties of bright X-ray emitting ellipticals, at least within their Re- 

Key words: galaxies: elliptical and lenticular, CD - galaxies: kinematics and dy- 
namics - galaxies: structure - galaxies: cooling flows - X-rays: galaxies - X-rays: ISM 



1 INTRODUCTION 

Since its discovery in the late 1970's (Fabbiano 1989) the 
presence of an extended and X-ray emitting interstellar 
medium (ISM) in early type galaxies raised the hope of 
having a new powerful and reliable tool to determine the 
dynamical properties of these galaxies, in particular their 
mass profile and total mass out to large radii (~ 10 optical 
effective radii -Re). In the simplest and commonly adopted 
approach, these determinations assume a spherical mass dis- 
tribution. Then, the most general expression for the momen- 
tum equation of the gas is: 

J^l^+V = ^g, (1) 
fxrUp p dr 

where T and p are the radial profiles of the ISM temperature 
and density, rrip is the proton mass, /i is the mean molecular 
weight and g = GM{r)/r^ is the total gravitational field. V 
includes all the non hydrostatic terms and, for simplicity, we 
do not consider other pressure terms such as the magnetic 
one (see however Sect. 4). Thus, in hydrostatic equilibrium, 
T(r) and p{r) - that in principle can be inferred from the 
X-ray surface brightness profile through deprojection, using 



also spectral information (e.g., Forman et al. 1985, Trinchieri 
et al. 1986) - are linked to the total mass profile by the well 
known formula (Fabricant, Lecar & Gorenstein 1980) 

M(.)^-.^f^ + ^). (2) 
Gfinip \ dlogr dlogr J 

Forman et al. (1985) made the first total mass determi- 
nations using this approach, for a sample of early type 
galaxies with Einstein data extending out to radii of ^ 50 
kpc. They obtained large mass values implying that dark 
matter should be present in a large quantity [M/Lb ^ 
25 — 100 Mq/Lb,©). However, uncertainties in the p and 
T profiles were so large that mass values could be even as 
low as those derived from optical measurements (Trinchieri, 
Fabbiano & Canizares 1986; Fabbiano 1989), that typically 
give M/Lb ~ 12Mq/Lb.q within Rc (Saglia, Berlin & Sti- 
aveUi 1992). 

Mass profiles can in fact be derived also from the ob- 
served stellar velocity dispersion profiles ap{R); these ex- 
tend out to ~ 1 — 2Rc in few best cases (e.g., Kronawitter 
et al. 2000). In this context, a common approach is to use 
the Jeans equations; for a spherical model described by a 
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phase-space distribution function / = f{E, L^), the relevant 
equation is: 



1 dp*(T^ ^ 2I3<tI 



dr 



= -9 



(3) 



(e.g., Binney & Tremaine 1987). Here p* is the stellar den- 
sity profile and /3(r) = 1 — erf (r)/2<T^(r), where (Tr is the 

radial stellar velocity dispersion and uf = + (t| is the to- 
tal tangential stellar velocity dispersion. /3(r) describes the 
radial trend of anisotropy in the stellar velocity dispersion 
tensor, which is radially anisotropic where < /? < 1, tan- 
gentially anisotropic where /3 < and isotropic where /3 = 0. 
For assumed visible and dark density profiles and a speci- 
fied radial trend of /3(r), eq. (3) is solved to derive tjf (r); 
then by projection one obtains <Jp{R) to be compared with 
the observed one. As well known, in absence of a physi- 
cally motivated choice for /3(r), this method suffers from a 
"dark matter-/? degeneracy" : a sufficient amount of orbital 
anisotropy can mimic the presence of a massive dark matter 
halo (e.g., Tonry 1983). Observationally, /3(r) can be derived 
from the deviation of the line profiles from a Gaussian (e.g., 
Gerhard et al. 1998), but this process requires high quality 
data and care in analysis. In fact, the use of the hot ISM to 
derive dynamical quantities w£is very welcomed also because 
it is exempt from this degeneracy problem^ . 

In hydrostatic equilibrium, the two conditions (1) (with 
V = 0) and (3), holding for the gas and the stars lying in the 
same total gravitational field g, allow us to determine M{r) 
and /3(r) uniquely, when coupled with the knowledge of the 
four observables T(r), p(r), p*(r) and ap{R). For example, 
by eliminating g from (1) and (3) one obtains a relation link- 
ing gas and stellar properties only (Fabricant et al. 1989): 



k dpT 1 dp,al 2/3a^ 



/urripp dr p* dr 



(4) 



Then, one can use the equation above to derive j3{r) using 
information coming from X-ray observations, as attempted 
recently by Mathews & Brighenti (2003a; hereafter MB03). 

The lack of reliable ISM density and temperature pro- 
files has been a severe limitation to the use of X-rays as a 
probe of the galactic potential and of stellar dyuaniical quan- 
tities. This was first overcome partially thanks to ROSAT 
observations. Currently, quite accurate T and p profiles are 
being derived from Chandra and XMM — Newton obser- 
vations (e.g., Sun et al. 2003). In addition, in a few cases 
Chandra has already shown signs of departure from the hy- 
drostatic equilibrium in the hot gas coronae (e.g., Jones et 
al. 2002). Therefore a timely question is how much reliable 
the X-ray estimate of dynamical quantities can be. 

In this paper we address the problem of the accuracy 
of the X-ray based total mass M(r) and /3(r) determination 
when the hydrostatic equilibrium hypothesis is assumed but 
it does not hold. Our analysis relies upon the use of the 
relevant hydrodynamical equations, while we do not discuss 
technical problems related with the extraction of physical 
quantities from the observations, such as the well known 



^ Note that in other approaches (e.g., Bertin et al. 1994), the 
anisotropy profile is not a function to be specified by independent 
arguments, but derives from the physical picture on which the 
model is based. 



ill-posed problem of Abel inversion for observed data (e.g., 
Tikhonov & Goncharsky 1987). 



2 ESTIMATES OF MASS AND ORBITAL 
ANISOTROPY WHEN V 7^ 

In the most general case (V ^ 0) the physically correct anar 

log of eq. (4) is 



k dpT _ 1_ dp^al _^ 2j3ijl 



lirupp dr 



p* dr 



(5) 



Instead, by erroneously assuming hydrostatic equilibrium, 
one uses an equation involving /S"^^* : 



k dpT _ 1 dp^al 2/3"Vf 



liTUpP dr p* dr 



(6) 



Therefore the difference between the estimated and the true 
anisotropy profile is given by 



(7) 



Following the same reasoning, the equation that should be 
used to derive the total mass profile M{r) is eq. (1), while 
using it with V = 0, one derives an estimated mass M"'* 
given by: 



M + -V. 



(8) 



Note that, although M*"* and derive from two indepen- 
dent procedures, the deviation in the mass estimate is linked 
to that in the /3 estimate by a remarkably simple relation: 



M-'(r) _ ^ 2ra^r{r)[rHr)-m] 
M{r) GM{r) 



(9) 



Thus, for example, M""* < M wherever (3""^ > 13. 

In general, the modulating factor 2r(j^/GM in eq. (9) 
depends on the exact slope of the stellar density profile, on 
the presence of dark matter and on f3{r) [via <Jr{r)]. The ra- 
dial behavior of this factor for a few representative cases is 
given in Fig. 1. Within ~ Q.2Re the factor is very sensitive 
to the shape of the density profile, as shown by the cases 
of the de Vaucouleurs (1948), .Jafi'e (1983) and Hernquist 
(1990) isotropic models (respectively the solid, long-d£ished 
and dotted lines). In addition. Fig. 1 also shows that this de- 
pendence is not trivially linked to the radial behavior of the 
p, profile, since p*(r) oc r""'^'*, r^^ and r^^ respectively 
for the de Vaucouleurs (Ciotti 1991), Jaffe and Hernquist 
models. The inner shape of the factor can be derived, how- 
ever, from the asymptotic analysis discussed by Bertin et al. 
(2002). 

The short-dashed line in Fig. 1 refers to an isotropic 
stellar de Vaucouleurs model coupled to a dark matter halo 
that is described by an Hernquist profile and is 5 times as 
massive as the stellar component; the mass enclosed within 
_Rc is the same for the two components (which implies a core 
radius for the Hernquist model of 2.47-Re). Compared to the 
corresponding pure stellar model (solid line) there is a small 
increase in the modulating factor. 
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When /3 > the factor is larger than in the isotropic 
case; for example, when /3 is radially constant^, as adopted 
in some modelings (see, e.g., Sect. 3), and the density dis- 
tribution is p* oc r~'^, the factor is 1/(7 — 1 — 0.5/3). In Fig. 
1 the dot-dashed line refers to a dc Vaucoulours profile with 
/9(r) radially varying according to the Osipkov-Merritt pa- 
rameterization (Osipkov 1979, Merritt 1985). In the central 
region, the effect of such a f3{r) = /{r^ + r^) (where Va 
is the anisotropy radius) is a small increase of the modu- 
lating factor with respect to the isotropic case (solid line). 
In the external regions an idea of the full range of the fac- 
tor values derives by comparing the cases of isotropic and 
anisotropic orbital distributions for a dc Vaucouleurs profile 
(solid and dot-dashed lines respectively). In fact, for this 
profile, the choice of ra ^ Re adopted for Fig. 1 gives the 
maximum allowed level of anisotropy, based on stability ar- 
guments (Ciotti & Lanzoni 1997; see also Sect. 3). 

In summary, outside the effective radius the modulating 
factor is everywhere < 0.7; for the "maximally" anisotropic 
case, it is everywhere < 1. Within 7?e it is again everywhere 
<1, except for the Hernquist model. Therefore for most mod- 
els "errors" on M (i.e., M''"'' /M - 1) are < than the /? de- 
viations (i.e., Z?""* - 13). 

In the next sections we estimate the difference /S*^** — /3 in 
eq. (7), first qualitatively (Sect. 2.1) and then quantitatively 
by using numerical models of hot gas flows for X-ray bright 
galaxies (Sect. 2.2). 

2.1 Qualitative analysis 

We assume that the hot ISM filling the galactic potential 
well has velocity, density and temperature u, p and T. Usu- 
ally, in the case of galactic flows, a mass source rate ap* (r) 
is present in the V term, due to the stellar mass losses (e.g., 
Ciotti et al. 1991). For an inviscid fluid, the V term is then 
written as 

V=^+l^ + ^. (10) 
at 2 dr p 

We consider first eq. (7) in the simplest case, namely the 
stationary case with a = 0; then 



oest :_ ^ 



r dv? 
4^ 



(11) 



This shows that when the velocity decreases with increasing 

radius in magnitude, i.e., when dv? /dr < 0, then /9'^''*(r) > 
P{r). This is the case of the central regions in cooling fiows. 
If, instead, the velocity increases in modulus with increasing 
radius, then /3(r) is underestimated. This could be the case 
of galactic winds. Note that fairly large differences fi""^* — (3 
are expected when v? varies by ~ over a characteristic 
length of the order of r. 

When a > 0, these results arc reinforced in two cases: 
when the gas is inflowing {u < 0) with a velocity de- 
creasing with radius {du'^/dr < 0), then even more radi- 
ally anisotropic orbits arc deduced; and when the gas is 
outflowing {u > 0) with a velocity increasing with radius 
{du^/dr > 0), then even more tangentially anisotropic or- 
bits are deduced. In the other cases the ap*/p term may 

^ Any consideration on the physical consistency of this assump- 
tion is left apart here; see, e.g., Camm 1952. 



have a competing effect on the du/dr term that depends on 
the magnitude of the source function a. 

All the results above can be converted straightforwardly 
into results concerning meiss estimates by using eq. (9) . 

The X-ray brightest galaxies, for which bettor quality 
data are available and dynamical studies can be attempted, 
contain large quantities of hot gas and are supposed to host 
inflows for which du^ /dr < 0. In their inner regions (well 
within ^ Re) the gas can reach high velocities, up to the 
local sound speed, while u ^ outside Re (e.g., Ciotti 
et al. 1991). Therefore, it is expected that more radially 
anisotropic orbits than true are deduced from X-ray obser- 
vations of the central regions; /S*"' derived outside Re should 
be reliable. More quantitative estimates follow below. 



2.2 Quantitative results 

Here we study the effect of the full V term given in eq. (10). 

For this investigation, two gas dynamical models are used, 
resulting from the numerical simulations of time evolving 
galactic gas flows by Pellegrini & Ciotti (1998) and represen- 
tative of hot gas rich galaxies. In principle, one should use 
a numerically computed right hand side of cq. (10). How- 
ever, in order to reduce the numerical noise due to spatial 
and temporal derivatives over a finite grid, the V term is 
obtained from eq. (1), so that: 



g + 



k dpT 
p.mpp dr 



Similarly: 



dpT 



M 



gpnipp dr 



(12) 



(13) 



In fact, for the galaxy models we know g and (5. As rep- 
resentative cases we consider 1) a bright inflowing galaxy 
{Lb = 7 X 10^° L©, Re = 9.0 kpc, central a = 287 km 
s~^ and a dark halo nine times as massive as the stellar 
component); 2) a less massive galaxy in which the inflow is 
confined within 0.1 Re {Lb = 2.2 x 10^° L©, Re = 4.4 kpc, 
central a = 210 km s^^ and a dark halo twice as massive 
as the stellar component). The underlying galaxy model is 
isotropic (/3 = 0) . 

Fig. 2 shows in the upper panel the radial trend of the 
gas velocity at the end of the simulations (corresponding to 
an elapsed time of 15 Gyrs). The panel below shows whether 
any term in the right hand side of eq. (10) gives a dominant 
contribution to V. Since the gas flows have reached a quasi 
steady state rRct 15 Gyrs, the du/dt term is negligible, and 
only the radial trend of the advective/source term ratio is 
plotted. For the largely infiowing galaxy the advective term 
is dominating in the central region and becomes comparable 
to the source term aX ^ Re, while for the partial wind model 
the two terms keep comparable. For both models the two 
terms produce a deviation in the same sense, since they are 
both negative in the inflowing region and positive in the 
outflowing region. 

In the lowest two panels. Fig. 2 shows what P""* and 
M'^"^ would be derived from X-ray observations. Overall, the 
trends of /J*^^* agree with the qualitative analysis of Sect. 2.1 
and, as expected, P"^"* is > for the central gas inflowing re- 
gion. For the galaxy hosting a global inflow, decidedly more 
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Figure 1. The modulating factor of cq. (9) for the isotropic Hern- 
quist (dotted line), JafTe (long-dashed line) and de Vaucouleurs 
(solid line) density profiles. The dot-dashed line represents a max- 
imally anisotropic dc Vaucouleurs model (see Sect. 2). The short- 
daslicd line is an isotropic stellar de Vaucouleurs model plus a 
Hernquist dark matter halo five times as massive as the stellar 
component; the mass enclosed within Rg is the same for the two 
components (see Sect. 2). 



radially anisotropic orbits than true would be deduced by ne- 
glecting the gas velocity; correspondingly, M^°* < M would 
be derived for the central galactic region [eq. (9) and Fig. 
2]. Outside the inflow region, P''^* is small and M''^* ~ M, 
consistently with the predictions for low gas velocities. 



3 THE CASE OF NGC4472 

The X-ray bright Virgo galaxy NGC4472 has been cho- 
sen since the earliest times for a detailed investigation of 
its dynamical properties using X-ray data (Thomas 1986; 

Trinchiori ot al. 1986; Locwonstcin 1992; Berlin, PignatoUi & 
Saglia 1993; Bcrtin & Toniazzo 1995; Irwin & Sarazin 1996). 
The study based on the most recent and extensive data set is 
that of MB03 (see also Mathews & Brighenti 2003b). MB03 
employed eq. (4) to relate /3(r) and T(r). They used p(r) de- 
rived from Abel-inverted Chandra, ROSAT and Einstein 
surface brightness profiles, T{r) derived from deprojectiori of 
ROSAT and Chandra data, the observed ap{R) from Ben- 
der et al. (1994) and Pried & Illingworth (1994) and p.(r) 
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Figure 2. For a global inflow (in a galaxy with Lb = 7.1 x 
IO^'^Lq, solid line) and a partial wind (in a galaxy with Lb = 
2.2 X IO^^Lq, dashed line; more details on the model galaxies are 
given in Sect. 2.2) from top to bottom the panels show the radial 

trends of: the gas flow velocity, the ratio between the advectivc 
and the source term in cq. (10), /S*^**, and M*^^* normalized to 
the true total mass. In the lowest panel the dotted lines give the 
stellar mass profiles of the two model galaxies (again normalized 
to the total mass; see Sect. 2.2 for more details). 



resulting from fitting the observed brightness profile with a 
de Vaucouleurs law (but in the central regions the density 
profile was assumed to steepen as r^" ''). 

Adopting a radially constant /3'^°* for the modeling, the 
value [3'^"'^ = 0.71 ± 0.15 was derived, while the optically de- 
termined /3{r) within R^, is everywhere <0.34 (Kronawitter 
et al. 2000). That is, (3""* is about twice as large as the max- 
imum (3 found from optical data. Conversely, assuming the 
optically determined 0{r) of Kronawitter ct al. (2000), T{r) 
higher than observed is obtained. Commenting this result, 
MB03 state that "This strong preference for radial stellar 
orbits must be understood in terms of the formation history 
of massive elliptical galaxies. Conversely, if the smaller, op- 
tically determined (3 is indeed correct, we are led to the im- 
portant conclusion that the temperature profile T{r) of the 
hot ISM in NGC 44"^^ must differ from that indicated by 



On the use of X-rays to determine dynamical properties of Elliptical galaxies 5 



X-ray observations, or that the hot gas is not in hydrostatic 
equilibrium." 

Here we argue that, if problems with the several steps 
involved in the modeling of the data can be excluded, the 
last possibility of lack of hydrostatic equilibrium is strongly 
favoured. In fact a simple stability argument rules out as 
unphysical such a high X-ray-derived /3^"* value. In case of 
a constant anisotropy parameter, the radial orbit instabil- 
ity indicator (Polyachenko & Shukhman 1981, Pridman & 
Polyachenko 1984) turns out to be 



2Kr 



1 



1 - 



3.45, 



(14) 



where Kr and Ki) + K^p are respectively the total radial and 
tangential kinetic energies. A galaxy is unstable due to radial 
orbit instability when this indicator is >1.7±0.25; therefore 
the value calculated in eq. (13) is well outside the stability 
regime. Note that even when inserting the maximum value 
of /3(r) = 0.34 derived by Kronawitter et al. (2000) in eq. 
(13), the indicator is « 1.5, safely in the stability regime. 
The criterion above has been proved to be obeyed essen- 
tially by all galaxy models to which it has been applied, 
even including massive dark halos (e.g., Stiavelli & Sparke 
1991; Nipoti, Londrillo & Ciotti 2002; Bertin & Trenti 2003; 
see also Saha 1991). 

Additional considerations axgue in favour of the hy- 
pothesis of lack of hydrostatic equilibrium to explain the 
inconsistencies between X-ray and optically-derived orbital 
anisotropy. In fact, there is also a deviation between the X- 
ray and the optically-derived mass for NGC4472, and the 
two discrepancies (concerning ,5 and M) are related as pre- 
dicted by eq. (9), as shown in the following. Having dis- 
carded the X-ray derived fl'^"* as unphysical, we assume that 
the optical estimate of (i is the correct one. Then, eq. (9) 
predicts the X-ray-derived M*^** to underestimate the true 
mass M. More specifically, the expected value of M"^"^ /M 
can be derived from the results for the de Vaucouleurs mod- 
els iu Fig. 1. For a range of values of the modulating fac- 
tor within Jie of ~ 0.5 - 1 and for - P = 0.4, then 
M^'^/M = 0.8 - 0.6 is given by eq. (9). What M"* derives 
from the X-ray data? The best fit total mass profile M*^"* (r) 
determined from eq. (2), with the same T{r) and p(r) used 
by MB03, is a one-component de Vaucouleurs mass profile 
out to r = Re, with a total constant mass-to-light ratio of 
W'^/Lb = 7 (Mathews & Brighenti 2003b). M/Lb esti- 
mates were derived also from optical studies based on high 
quality photometrical and kinematical data extending out 
to ~ 0.5iie. The M/Lb determined from one-component 
constant mass-to-light ratio axisymmetric dynamical mod- 
els is 9.2, for the central galactic region, when rescaled to 
the distance adopted by MB03 (van dcr Marel 1991). An- 
other analysis used spherical self-consistent two-component 
models (Saglia et al. 1992). The models selected at a 3 sigma 
confidence level have M/Lb = 10.1 — 14.6 within _Ro, when 
rescaled to the MB03 distance'*. These optically determined 
M values, if taken as true values, imply that M*^"* underesti- 
mates the true mass M. More precisely, M^"*- /M = 0.76 and 



^ These selected models also have a velocity dispersion profile 
consistent with the data for the globular cluster system out to 
~ 2.5ile. 



0.69 — 0.48 when adopting respectively the van der Marel's 
and the Saglia ct al.'s results, in remarkable agreement with 
the predictions above {M""^ jM = 0.8 - 0.6). Therefore, the 
hypotheses of i) correct optical estimates of /3 and M and 
ii) departure from hydrostatic equilibrium of the hot ISM, 
give a consistent explanation for the discrepancies of the 
X-ray-derived (3^^* and M""* values. In particular, these hy- 
potheses account for the fact that is an overestimate 
and M''"* is an underestimate and also for the size of the 
discrepancies with the true (3 and M. 



4 DISCUSSION AND CONCLUSIONS 

We have investigated the efi'ect that deviations from hydro- 
static equilibrium have on the X-ray estimates of the to- 
tal mass M and the orbital anisotropy parameter f3 via the 
"standard" eqs. (2) and (4). We have found that: 

• in the most general CEise, the deviations of /J*^^* and 
M"^* from the true values are linked by a remarkably sim- 
ple relation. In particular, the mass is underestimated if /3 
is overestimated. Also, for the most common density and 
orbital anisotropy distributions, with or without a massive 
dark halo, the /J*^"* deviations are larger than the "errors" 
on M at all radii >0.2i?e. 

• in the stationary case without mass/momentum 
sources/sinks, more radially arnsotropic orbits than true are 
deduced for the central regions of gas inflowing galaxies. In- 
stead, /3(r) is underestimated for outflowing regions where 
the velocity increases with radius, as in galactic winds. When 
mass sources are present, the results above are reinforced in 
the inflow and in the wind cases. In the other cases the mass 
source term may have a competing efi'ect. 

• the full time dependent problem is investigated with 
the aid of numerical simulations reproducing gas rich galax- 
ies. The results are consistent with the simple qualitative 
predictions. /J*^^* decidedly overestimates the true value (and 
M"^* underestimates it) if derived from the central region 
where gas is inflowing with a large velocity, while outside it 

~ f3 and M"* - M. 

• and derived for the most thoroughly inves- 
tigated bright eUiptical galaxy (NGC4472), for which the 
best X-ray data are available, are considered as a repre- 
sentative application of our results. We show that the high 
(3""* recently derived is unphysical. We also note a discrep- 
ancy in the X-ray and optically derived masses within Re. 
Assuming as true /3 and M the optically derived values, 
the /3-overestimate implies (when applying our results for 
the non-hydrostatic equilibrium) an M-underestimate in re- 
markable agreement with M*^'* — M. Of course, one cannot 
exclude that a different modeling of the X-ray and optical 
data would be able to produce a /3(r) profile consistent with 
the observed one, in the hydrostatic equilibrium hypothesis. 
However, this would reveal a degeneracy in the procedure 
and undermine the possibility to uniquely derive dynamical 
information from the X-ray data. 

• note that eqs. (7-9) describe the deviations of /J*^"* and 
M"^* from the true values even when sources of pressure ad- 
ditional to the thermal one are effective (such as magnetic 
and radiation pressure). In fact, if these are accounted for 
explicitely in eqs. (2) and (4), then eqs. (7-9) remain unal- 
tered. If these are not accounted for, again eq. (9) remains 
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unaltered and eqs.(7-8) hold provided that the additional 
pressure terms are considered as included in the V term. 
The deviations of /S*^"* and M*^*' produced by these kinds 
of pressure may be confined to the central galactic regions, 
where magnetic fields and AGN-related phenomena more 
likely are concentrated. 

Additional explanations can be figured out for the devi- 
ations of X-ray-derived (3""* and M"^^* from the true values. 
These include deviations from spherical symmetry, projec- 
tion/deprojection eff'ccts not adequately taken into account, 
presence of a multiphase ISM at each radius, amplification 
of uncertainties in the process of extracting T{r), p{r), p*{r) 
and Orir) from the data. However, it is yet to be investigated 
the possibility for these explanations to account at the same 
time for deviations in /3 and M, and by precisely the ob- 
served amount, as demonstrated here for the hypothesis of 
lack of hydrostatic equilibrium in the case of NGC4472. 

Our analysis casts doubts on the possibility to use X-ray 
information and the hydrostatic equilibrium assumption as 
robust tools to determine accurately the dynamical proper- 
ties of X-ray emitting galaxies, at least within Re- Another 
warning to take with care X-ray mass estimates had already 
been made when it appeared that mass models constrained 
by accurate and extended stellar dynamical data, and in- 
cluding different amounts of dark matter, all lead to simi- 
lar temperature and X-ray brightness profiles (Bertin et al. 
1993), within the uncertainties in the X-ray measurements. 
In this sense, the situation might be improving with the new 
Chandra and XMM data. 

In conclusion, instead of obtaining robust dynamical es- 
timates from hot gas-based studies, it should perhaps be ex- 
pected that stellar dynamical studies (possibly coupled with 
gravitational lensing studies) help understand the dynamical 
status of the hot gas. 
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